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Abstract 

The minimum mass that a virialized gas cloud must have in order to be 
able to cool in a Hubble time is computed, using a detailed treatment of the 
chemistry of molecular hydrogen. With a simple model for halo profiles, we 
reduce the problem to that of numerically integrating a system of chemical 
equations. The results agree well with numerically expensive 3D simulations, 
and our approach has the advantage of rapidly being able to explore large 
regions of parameter space. The minimum baryonic mass M& is found to be 
strongly redshift dependent, dropping from 10 6 Mq at z ~ 15to5xlO 3 M at 
z ~ 100 as molecular cooling becomes effective. For z S> 100, rises again, 
as CMB photons inhibit ^-formation through the H~ channel. Finally, for 
z ^> 200, the /^-channel f° r -^-formation becomes effective, driving Mb 
down towards Mf, ~ 10 3 Mq. With a standard CDM power spectrum with 
cjg = 0.7, this implies that a fraction 10 -3 of all baryons may have formed 
luminous objects by z = 30, which could be sufficient to reheat the universe. 
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1 INTRODUCTION 



1.1 When did the universe reheat? Observational constraints 

It is now widely accepted that the universe underwent a reheating phase 
at some point after the standard recombination epoch at redshift z ~ 10 3 . 
However, the question of when this happened remains open. The absence of 
a Gunn- Peterson trough in the spectra of high-redshift quasars has provided 
strong evidence for the reheating occurring at a redshift z > 5, since it 
indicates that the intergalactic medium (IGM) was highly ionized at lower 
redshifts (Gunn and Peterson 1965, Steidel and Sargent 1987, Webb et al. 
1992). The smallest baryonic objects to go non- linear in a standard cold 
dark matter (CDM) model are expected to reionize the IGM at a redshift 
somewhere in the range 10 < z < 100 (Bond k Szalay 1983; Couchman 
1985; Couchman k Rees 1986; Fukugita k Kawasaki 1991; Tegmark, Silk 
k Blanchard 1994; Tegmark k Silk 1995; Liddle k Lyth 1995). In recent 
models with baryonic dark matter, reheating and reionization is predicted 
to occur at an even higher redshift, typically in the range 100 < z < 1000 
(Peebles 1987, Gnedin k Ostriker 1992, Cen et al. 1993). 

A reheating epoch would have at least two interesting classes of effects 
that may be measurable today: effects on subsequent structure formation 
and effects on the cosmic microwave background radiation (CMB). Subse- 
quent structure formation would be affected in at least two ways: 

1. The heating of the IGM up to a higher adiabat would raise the Jeans 
mass, thus suppressing the formation of small objects. For instance, 
an IGM temperature of 10 5 K at a redshift of a few would suppress 
the formation of galaxies of mass below 10 10 M Q , thus alleviating the 
ubiquitous problem of theories overpredicting the abundance of faint 
galaxies (e.g. Blanchard et al. 1992, Kauffman et al. 1993; Cole et al. 
1994). 

2. If the objects that reheat the IGM also enrich it with heavy elements, 
the ability of gas to cool would be greatly enhanced in the temperature 
range 10 4 K < T < 10 7 K, presumably facilitating future structure 
formation. 

The CMB would be affected in at least three ways: 

1. Hot ionized IGM would cause spectral distortions which might violate 
the stringent limits on the the Compton y-parameter (Mather et al. 
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1994) . This is a problem mainly for BDM models (Tegmark & Silk 

1995) . 

2. Spatial fluctuations on angular scales below a few degrees may be 
suppressed, while fluctuations on larger scales would remain fairly un- 
affected. Therefore a comparison of the results of current and future 
degree scale experiments with those of COBE (Smoot et al. 1992) con- 
strains the reionization epoch. 

3. New spatial fluctuations will be generated on smaller angular scales, 
through the so called Vishniac effect (Vishniac 1988, Hu et al. 1994). 
The current upper limit on CMB fluctuations on the 1 arcminute scale 
(Subrahmanyan et al. 1993) places constraints on some reheating sce- 
narios. 

In other words, with the recent surge in CMB experiments and the consider- 
able numerical, theoretical and observational results on structure formation, 
the thermal history of the universe is now coming within reach of our exper- 
imental probes. In view of this, it is very timely to theoretically investigate 
the nature of the reheating epoch in greater detail, and investigate the prop- 
erties of the objects that caused it. In this paper we will focus on two of the 
most basic attributes of these first objects: their mass and their formation 
redshift. Hence the goal is to derive the mass-redshift distribution of the 
very first objects that might be able to reheat the universe, and thus set the 
stage for all subsequent cosmological events. 

1.2 What does theory predict? 

In both CDM and BDM models of structure formation, the first objects 
predicted to go nonlinear are the smallest ones. The crucial question is if 
cooling will allow the baryonic clouds to dissipate their kinetic energy and 
collapse more than the dark matter, to eventually become self-gravitating 
and form an interesting object like a galaxy, a V.M.O. or a black hole (see 
e.g. Binney 1977, Rees and Ostriker 1977; Silk 1977; White and Rees 1978; 
Araujo & Opher 1988, 1989, 1991). For low mass objects, the smaller they 
are, the less efficiently they dissipate energy and cool. Thus a detailed 
treatment of gas-dynamical processes will predict a characteristic mass scale 
M c such that objects with M > M c can cool rapidly, whereas smaller lumps 
will merely remain pressure supported and not form anything luminous. In 
other words, M c is the mass scale of the first luminous objects. 
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Fortunately, making a theoretical estimate of M c is much simpler than 
the corresponding problem for present-day structure formation. Today there 
are large uncertainties in both the metal abundance of the intergalactic 
medium (IGM), which affects cooling rates, and in the UV background, 
which affects ionization rates and molecular chemistry. Before the first struc- 
tures formed, there was by definition neither metals nor UV background. 

The problem has recently been treated realistically using a multi-fluid 
3D cosmological hydrodynamics code which evolves not only the dark, and 
baryonic matter, but also tracks the non-equilibrium chemistry of 9 species, 
including hydrogen molecules (Abel 1995; Anninos et al. 1996; Norman et 
al. 1996). The main obstacle to this program is computational expense, 
because of the large dynamical range involved. As a complement to such 
heavy computations, it is thus worthwhile to attack the problem with vari- 
ous approximate techniques that are fast enough to run many times, thereby 
exploring all of parameter space and finding out which parameter choices 
and initial conditions merit more detailed numerical studies. This is the 
purpose of the present paper. One such approximate method is that of 
Haiman, Thoul & Loeb (1996, hereafter HTL96), which numerically fol- 
lows the growth if an isolated density peak that is spherically symmetric. 
Although the first structures to collapse in CDM are typically sheet-like 
rather than spherically symmetric, this model nonetheless illustrates which 
physical processes are likely to be the most important in the full 3D case. 
Since the approach of HTL96 involves numerically integrating a partial dif- 
ferential equation (separately tracking a large number of spherical shells), 
it is still fairly time-consuming, and results are presented for only 24 points 
in the M — z plane (see Figure ||). In this paper, we use a still simpler ap- 
proach, involving nothing but ordinary differential equations, which turns 
out to reproduce the results of HTL96 quite well. The resulting code is 
so fast that we can run it thousands of times, thereby finding the curve in 
Figure ^ that delimit collapsing objects from non-collapsing ones, and study 
how this curve depends on cosmological parameters such as fi, fi& and h. 

For the various CDM-based scenarios, the first interesting objects will 
turn out to be rare peaks in the Gaussian random field of mass between 10 4 
and 10' Mq, at redshifts in the range 20 < z < 100. At these redshifts, the 
initial IGM temperature is considerably lower than the virial temperatures in 
question, so the baryons will initially collapse together with the dark matter. 
These first objects will have virial temperatures between a few hundred and 
a few thousand degrees, which means that the main coolant will be molecular 
hydrogen. (Line cooling by hydrogen and helium is negligible for T <C 10 4 K, 
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and lithium hydride and other less abundant molecules become dominant 
only when T <C 500K.) In Section 2, accurate expressions for H2 cooling are 
presented, and it is shown that the pre-collapse H2 abundance is typically 
too low for the clouds to cool significantly in a Hubble time. The fate of a 
virialized lump thus depends crucially on its ability to rapidly produce more 
H2, which is the topic of Section 3. Our simple model for the evolution of 
density and temperature is presented in Section 4, and the numerical results 
are described in Section 5. The results and their cosmological implications 
are discussed in Section 6. 

2 COOLING BY MOLECULAR HYDROGEN 

How much molecular hydrogen is needed for a gas cloud to be able to cool in 
a Hubble time? This question will be answered in the present section. The 
atomic physics of molecular hydrogen cooling has been studied extensively 
by many authors, e.g. Lepp & Shull (1983). An excellent review of what will 
be needed here is given by Hollenbach and McKee 1979 (hereafter HM), who 
also provide a number of useful analytical fits to various numerical results. 

When an H2 molecule gets rotationally or vibrationally excited through 
a collision with an H atom or another H2 molecule, there are two competing 
channels through which the ensuing de-excitation can occur. Either the de- 
excitation is radiative, which amounts to cooling, or it is collisional, in which 
case there is no net energy loss from the gas. When the density n is very low, 
the former channel dominates. In this case, the hydrogen molecules spend 
most of their time in the ground state or in the J = 1 rotational state (whose 
radiative decay to the ground state is forbidden since H2 has no dipole 
moment), and collisional excitations are for all practical purposes instantly 
followed by a radiative decay. Thus in the low density limit, the energy loss 
per unit volume is proportional to n 2 . When the density n is very high, 
on the other hand, collisions dominate. Thus to a good approximation, the 
distribution of molecules in various states is the the Boltzmann distribution 
of LTE, local thermal equilibrium, and the energy loss per unit volume is 
only linear in n. The border between "high" and "low" density is roughly the 
function n cr defined below. It is temperature-dependent, but lies between 
10 3 and 10 4 cm~ 3 for our regime of interest, 10 2 K < T < 10 3 K. A just 
virialized gas cloud has an overdensity of about 18-7T 2 « 178, i.e., a hydrogen 
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density 



n«23cm- 3 (^U 1 3 00i (i) 



so during the early stages of collapse, we are well into the low-density regime 
for our parameter range of interest. (Here and throughout this paper we 
assume a Helium abundance of 24% by mass.) 

Since the fraction / of hydrogen in molecular form will be quite low in 
our application (typically below 10~ 3 ), we can neglect H2 — H2 collisions 
and the formulas of HM reduce to the following: The cooling rate is 



t (He) 

L sa Lr - (2) 

1 + n cr /n 



where the critical density is 



L Qte) 

which depends only on temperature, not on n. Here the cooling rate in LTE 
is 

L (lte) w I ( ( 9 ^ Xl n °"^3' 76 > \ e -(0-13/T 3 )3 + 3 x l -24 e -0.51/T 3 | ergcm 3 s -l ; 



n 



1 + 0.12T3 2 - 1 



(4) 

whereas the cooling rate in the low-density limit is 

4-0) w ^ 2{E2 _ Eo)e -(E 2 -E )/ k T + ?_ Tj{E3 _ Ei)e -(E s - El )/kT_ (5) 

Here T 3 = T/1000K, Ej = J(J + l)E 1 /2, where Ek/k ~ 171K. Thus (E 2 - 
E )/k = (3/5) (J5 3 - Ei)/k = 3E 1 /k w 512K. 72 and 73 are the collisional 
de-excitation rates from the J = 2 and J = 3 rotational levels. The rates for 
collisional quadrupole de-excitation J — > J — 2 due to impact of hydrogen 
atoms are well fit by (HM) 



fl()- n T 1/2 \ ( 

7j(T) = l i + 6or- 4 + 10 " 12T3 J (°- 33 + a9ex P 



J -3.5 
0.9 



cm 3 s 1 . 



(6) 

Equation (||) assumes an ortho-H2 to par a- H2 ratio of 3:1. The first term 
gives the cooling contribution of para-H 2 and the second that of ortho-H2. 
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Abel et al. (1996b) show that for H2 formation by the gas phase reactions 
discussed in the following section, the interconversion mechanism, 

R 2 (ortho) + H+ -> R 2 (para) + H+ (7) 

will be fast enough to convert all ortho-H2 to para-H2 . Hence the appropriate 
cooling rate is given by the first term in equation (||) multiplied by four, i.e., 
equation (||) is replaced by simply 

L^°) « 5 l2 (E 2 - E )e-^- E ^/ kT . (8) 

Defining the cooling timescale t coo \ = T/T, we thus obtain 

r*** 48200 years (l + ^^(M)" 1 , (9) 

where n\ = n/lcm" 3 . Let us define the Hubble timescale at a redshift z 
as the age of the universe at that redshift. Then for = 1, 

r h ps 6.5 x 10 6 years h~ l z^ 2 . (10) 

Since the primordial gas clouds in which we are interested have just virial- 
ized, the Hubble timescale is of the same order as the gravitational timescale 
T g = (pG)~ l l 2 , the timescale on which collapse would proceed if the tem- 
perature where lowered and the clouds lost their pressure support. Thus 
the future of a newly formed gas cloud is crucially dependent on the ratio 
Tcooi/Th (Rees & Ostriker 1977). If r coo \ -C Th, the gas cloud will rapidly cool 
and begin a nearly free-fall collapse, whereas if t coo \ ^> t/j, the cloud will 
remain pressure supported and fairly stationary until much lower redshifts. 
Tcool = Th for 

f « 0.00016 ( zj 2 1 + i^-j ) e 512K / T . (11) 
J V0.03/ 100 y 60 + T 3 4 y v ' 

This critical H 2 fraction is plotted in Figure [j], as a function of temperature. 
It is seen that the H 2 fraction required exceeds typical initial abundances 
(~ 10 -6 ) for all redshifts z < 200 when T < 10 4 K. Thus our low-mass, 
high-redshift clouds can cool and collapse only if additional H 2 is produced 
(unless T ^ 10 4 K, in which case hydrogen line cooling will be effective). 

In the following section, we will compute how much additional H 2 will be 
produced, and discuss the conditions for when sufficient cooling will indeed 
occur. 
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Reaction 


Kate fc [cm /sj 


I \ (M CI (M let 1 


H+ + e" i-> H + hv 
H + e~ i-> H + /iz/ 

r+HMH 2 + e - 

H + |Hh + /V 
H^ + HmH 2 + H+ 


h Rs 1.88 x 10- 1U T- U (34 
fc 2 « 1.83 x 10- 18 T - 88 
fc 3 « 1.3 x 10~ 9 
k 5 w 1.85 x 10- 23 T 18 
A; 6 « 6.4 x 10" 10 


Hutchins (1976) 
Hutchins (1976) 
Hirasawa (1969) 
Shapiro & Kang (1987) 
Karpas et al. (1979) 


Reaction 


Rate fc [1/s] 


Reference 


H" + /iz/ i-> H + e~ 
H^ + hv i-> H + 

e~ + i— > e~ + 


fc 4 « 0.114T^ 13 e- 865(J /^ 
k 7 « 6.36 x io 5 e - 7160 °/ T ^ 
£ 8 » 4.91 x 10" 22 T4 


Appendix A 
Appendix A 



Table 1: Reaction rates used (all temperatures in Kelvin) 



3 PRODUCTION OF MOLECULAR HYDROGEN 

How much molecular hydrogen will be produced in a Hubble time? In 
hydrogen of density n = n[H] + n[H + ] + 2n[i?2] a t temperature T < 10 3 K, 
the ionization fraction x = n[H + ]/n and the molecular fraction / = n[H2]/n 
evolve as 

x = —k\nx 2 , (12) 

/ = k m n(l -x- 2f)x. (13) 

Collisional ionization of H atoms as well as collisional dissociation of H2 
is completely negligible at such low temperatures. At the low densities in 
question, H2 is formed mainly via the reaction H + e~~ — > H~ + hv at the 
rate fo, after which one of the following two things happen to the H~ almost 
instantaneously: 

• Molecular hydrogen is produced through the reaction H + H~ — > H2 + 
e~, at the rate k^. 

• The H~ gets destroyed by a CMB photon, at the rate k±. 

Thus the effective rate of information is k2kz/\k?, + k±/{\ — x)n\. Since 
T 7 cx 1 + z, the exponential term in £4 effectively makes ^-production 
through the //"-channel impossible for z S> 200. A second, less effective 
channel for molecule formation is the slow reaction H + + H — ► + hv at 
the rate /C5, followed almost immediately by either + H —> H2 + H + at 
the rate k% or photodissociation at the rate kj, thus producing H2 at the 
net rate k^k^/lke + kr/(l — x)n\. This channel works up to higher redshifts, 
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but since fcs <C k 2 , it becomes important only for lumps with virialization 
redshifts z V i T S> 100. In our calculations further on, we use the exact rate, 
i.e., 



k n 



k-A 



k% + ki/ (1 — x)n 



k 2 + 



ke + fa /{I — x)n 



k 5 . (14) 



Although we integrate the above-mentioned chemical equations numerically 
in our analysis, a number of the features of the solutions can be readily 
understood from the following elementary observations. First note that 
equation (jl^) is independent of /, since the electrons act only as a catalyst 



in the reactions that produce H<i- Since the right-hand-side of equation (12) 
is not linear but quadratic in x, the residual ionization fraction decays much 
slower than exponentially. In the absence of cooling, T and n will remain 
roughly constant in the pressure-supported cloud, and the solution will be 

x(t) = 5°-— , (15) 

w 1 + x Q nfat ' v ; 

i.e., x — > only as 1/t, where fa is the recombination ratej^j Substituting 
this into equation (|l3|), we see that / — > 1 as t — > oo, i.e., all hydrogen would 
become molecular if we waited long enough. With parameters in our range 
of interest, however, / will remain much less than unity for many Hubble 
times. Thus taking 1 — x — 2/ ~ 1, equation flT3| ) has the solution 

f(t) = fo + ^Hl + xonht) (16) 
fa 

when k m is roughly constant (it will be roughly constant except at z ~ 300 
and z ~ 100, which is when the two radiative dissociation processes go from 
being dominant to negligible). Thus the time evolution separates into two 
distinct regimes: XQnfat <C 1 and x^nfat 1. In the first regime, the 
residual ionization remains roughly constant, and molecules get produced 
at a constant rate. In the second regime, electron depletion becomes a 
serious problem, and the molecular fraction grows only logarithmically with 
time. Since the factor \/{xonk) is simply the recombination timescale, we 
can rephrase this result as stating that the molecule fraction produced is 
f — fo = (k m /ki) ln(l + N rec ), where N rec is the number of recombination 
times elapsed. The transition occurs after about one recombination time 



2 To obtain better accuracy when z ~ 10 3 , we use the more complicated rate equations 
given in Peebles (1993, §6) in place of the rate k\ from Table 1 in our numerical runs. 
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_l -\2.13 -3173/(1+*) 



-1 



(N rec ~ 1), i.e., when 

/~/ c = ^~ 3 . 5x 10- 4 T 3 L52 [l + 7.4 x 10V A (1 + zf A ' 6 e 

(17) 

for z <C 300, a value that is independent of the initial ionization fraction 
Xq. The factor in square brackets corresponds to photodissociation of H~ , 
and can be ignored for z < 100. Figure |l] shows /(t/j) as a function of T 
for xo = 3 x 10" 4 , together with / C (T). As can be seen, we typically have 
f{ T h) > fc for z v i r S> 50, i.e., we are well into the electron depletion regime, 
which means that the final molecule abundance / is rather insensitive to the 
initial ionized fraction xo and approximately given by equation (fTij). 

Figure U also shows that the three solid dots almost line up horizontally. 
In other words, the molecular fraction in clouds that just barely manage to 
collapse (where the molecular hydrogen fraction produced within a Hubble 
time is just enough to make it cool in a Hubble time) is almost independent 
of the virialization redshift for 25 < z v { r < 100. Since the virial temperature 
of a collapsing cloud is determined only by its mass and its virialization 
redshift, this implies that any cloud with a molecular hydrogen fraction 
~ 5 x 10~ 4 is able to cool within a Hubble time. We can summarize this 
with the following useful rule of thumb: If the virial temperature is high 
enough to produce a molecular hydrogen fraction of order 5 x 10~ 4 , then 
the cloud will collapse. This explains the rather constant slope in figures || 
and E for 20 < z mr < 80. 



4 EVOLUTION OF DENSITY AND TEMPER- 
ATURE 

In this section, we describe our simple model for how the gas density and 
temperature evolve in an overdensity that grows, goes nonlinear and virial- 
izes. Section 4.1 refers mainly to the dark matter — the late stages of the 
density evolution of the baryons is discussed in 4.2 and 4.3. 



4.1 The density 

Early on, while z 3> > space is approximately flat and the Friedmann 
equation has the approximate solution 

a{t) cx t 2 / 3 (18) 
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T=<4^ +oiaS,3) - <2i » 



regardless of the values of and the cosmological constant Ao- If an = 1 
universe has a completely uniform density p except for a "top hat" overden- 
sity, a spherical region where the density is some constant p' > p, then this 
top hat region will gradually begin to expand slower than the rest of the 
universe, stop expanding and recollapse to a point. By Birkhoff's theorem, 
the radius of this region will evolve according to the Friedmann equation, 
but with some O > 1. It is well known that the overdensity 

6=--l (19) 
P 

evolves as 

« 9 (a — sina) 2 ., 3 9 , , n . 

2 (1 — cos a) 6 20 

where the parameter a, the "development angle" is related to the redshift 
through 

1 + z v i r fa — sin a \ 2 / 3 a 2 
T+z = V 2^ ) = (127)2/3 

Here z v i r is the redshift at which the top hat would collapse to a point. In 
reality, an overdense region would of course not collapse to a point (and 
form a black hole). Since it would not be perfectly spherically symmetric, 
collisionless dark matter particles would mostly miss each other as they 
whizzed past the central region and out again on the other side, eventually 
settling down in some (quasi-) equilibrium configuration known as the virial 
state. For baryons, gas-dynamical processes become important, and pressure 
eventually halts the collapse at some density p p as discussed in Section [D| 
below. Strictly speaking, virial states are not stable over extremely long 
periods of time, and their density is certainly not uniform. For a virialized 
lump, often referred to as a "halo", a typical density profile peaks around 
some constant value in its core and falls off like 1/r 2 over some range of 
radii. Nonetheless, halos are often said to have a "typical" density 

p viT ^l8TT 2 p (l + z mr ) 3 , (22) 

which is a useful rule of thumb. Thus in the top-hat collapse model, density 
in the perturbed region is assumed to evolve as in Figure ||: the density 
starts out decreasing almost as fast as the background density p, with 

5 oc (1 + z)- 1 
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early on, just as in linear theory, but gradually stops decreasing and increases 
radically as z approaches z v i r . It never increases past the virial value p v ir 
or the pressure-determined value p p , whichever is smaller, but stays at that 
density for all z < z v i r . The main motivation for the use of the Lagrangian 
code in HTL96 was to provide a more realistical modeling of the spatial 
structure of the halo. We use the simple top-hat approximation instead, for 
the following reasons: 

• It requires much less computer time. 

• It reproduces the results of HTL96 fairly well. 

• The spherical symmetry assumption of HTL96 is probably somewhat 
inaccurate anyway, since n-body simulations have demonstrated that 
the first collapsed structures tend to be sheetlike pancakes rather than 
spherically symmetric. 

In defense of the spherical symmetry assumption, very rare peaks in a ran- 
dom field (which might correspond to the very first objects) are typically 
almost spherically symmetric (Bardeen et al. 1986). More importantly, since 
the virial temperatures in our application are typically only slightly higher 
than the pre-collapse gas temperatures, none of our conclusions should be 
very sensitive to the actual way in which the cloud gets to its virial configu- 
ration, such as whether it first passes through an intermediate pancake-like 
configuration or not. 

Unfortunately, a cannot be eliminated from the equations that relate 5 
and z by using elementary functions. For this reason, we use the following 
fit to the density evolution p(z) = p [l + S(z)], which is accurate to about 
5% until z is within 10% of z v i r (Tegmark 1994), at which the density is 
assumed to start approaching the limiting value p v - u or p p anyway: 



and po is mean density of the universe today. We use this fit in our numer- 
ical analysis, but never let the density exceed the virial value, as shown in 




(23) 



where 




(24) 



Figure §. 
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4.2 The temperature 

The thermal evolution of the gas is dominated by the following processes: 

• Hydrogen line cooling (as given by equation (pq)) 

• Cooling by molecular hydrogen (as given by equation (^)) 

• Compton cooling (as given by equation (pq)) 



• Adiabatic cooling/heating (caused by the expansion/compression of 
the gas) 

Bremsstrahlung and helium line cooling are completely negligible at the low 
temperatures in which we are interested. The first three mechanisms simply 
couple the gas atoms to the radiation field, which means that they will 
cause cooling when the gas is hotter than the CMB and heating otherwise. 
In other words, none of these mechanisms can make the gas cooler than the 
CMB temperature, which at z = 100 is a few hundred K^| In the Compton 
case, this is reflected by the fact that the cooling rate is of the form 



dT\ 
~dt) 



= k s x(T 1 -T). (25) 

comp 



For line cooling, given by (Dalgarno & McCray 1972) 

A; « 7.5 x 10- 19 erg cmV V 118348K / T n 2 *(l - x), (26) 

the CMB temperature is completely irrelevant, since line cooling only be- 
comes important when T S> 10 3 K, i.e., when T ^> T~. For the molec- 
ular case, this is included by replacing A m (T) by the net cooling rate 
A m (T) -A m (T 7 ). 

The adiabatic contribution is given by the pdV work done as the gas 
expands or contracts. In the simple top-hat model of the previous section, 
the density of the lump remains almost uniform until close to the virialization 
redshift z v i r , so that the adiabatic cooling term is simply 



dT\ 2n m 

-77 ) =o~T, (27) 

at J adiab & n 



3 Assuming nucleosynthesis abundances, cooling by lithium hydride is negligible com- 
pared to _ff2-cooling unless T <C 100K (Puy el al. 1993, Puy & Signore 1995), so we can 
safely neglect lithium chemistry for our application. 
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where the baryon number density n oc p is given by equation (^). (The 
molecular abundances are so small that to a good approximation, we can 
treat the IGM as a 7 = 5/3 monoatomic ideal gas.) As z — > z v i r , equa- 
tion (|23| ) would imply that T — > 00, as the lump collapses to a point. Instead, 
the lump is assumed to settle into an approximately pressure-supported con- 
figuration, where a typical gas element will obtain the virial temperature 
T v i r . For an overdense lump of total (baryonic and dark) mass M that stops 
expanding, recollapses and virializes at redshift z v i r , this temperature T„j r , 
which corresponds to the gas particles having similar velocities as the dark 
matter particles, is approximately (Blanchard et al. 1992) 

4.3 The effect of gas pressure 

How high will the typical gas density be in this pressure-supported state? At 
redshifts 3> 100, the Compton coupling to the CMB via the small fraction 
(10 -5 — 10 -3 ) of the electrons that remain ionized is still so strong that the 
IGM temperature will be close to that of the CMB, 

T 7 « 273K (^p) • (29) 

As time progresses, the Compton coupling weakens, and the IGM begins 
to cool below the temperature T 7 , cooling adiabatically as (1 + z) 2 . Com- 



paring equation (28) and equation (^), we therefore see that as long as 
M ^> 10 4 Mq, the baryons in the ambient IGM will have a temperature con- 
siderably below T v i r , and begin to fall into a virial configuration together 
with the cold dark matter. However, the gas density can only rise by the 
large factor 187r 2 without problems with pressure support if T <C T v i r after 
the collapse. Since T oc n 2//3 during the adiabatic compression, this means 
that W6 must have T V i r > (18vr 2 ) 2 / 3 T ~ 32 T before the collapse to be 
able to ignore pressure, and this turns out to be a good approximation for 
the critical masses only when z v i r <C 100. Otherwise, the condition that 
T = T v i r after the collapse gives only a collapse factor of order (T^ r /Ti) 3 / 2 , 
where T% denotes the temperature of the uniform background medium at 
redshift z = z v i r . In other words, equation (p2|) is replaced by 



T ■ \ 3//2 

p p » po(l + zwrf ( -jr 1 ) ■ (30) 
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We can obtain a more rigorous estimate of the final density as follows (Loeb 
1996). Hydrostatic equilibrium after the collapse implies that gravity is 
balanced by pressure gradients, i.e., that the gravitational potential <f> and 
the pressure p are related by 



V</> 



--Vp. 

p 



(31) 



Integrating this equation along some curve from very far outside the lump 
(where 4> = by definition) to a typical point inside the lump, we thus 
obtain 

Vp 



(32) 



Since the gas has been compressed adiabatically during the collapse to this 
state, its pressure and density are related by 



^5/3 
Pi) 



(33) 



where p\ and p\ denote the pressure and density of the uniform background 
medium at redshift z = z v i r . Substituting this into equation (32), we obtain 



\ 2 / 5 



OX = — - — 



5pi 
2 pi 



p_\ 
Pi J 



2/5 



1 



(34) 



By the ideal gas law, pi/pi = kTx/m p , where m p is the molecular weight. 
Eliminating (p/pi) using equation ( |33[ ) and defining T v i r by 



1 

--m z 
2 p 



we thus find the final overdensity inside the lump to be 



(1 + 5) 



P_ 

Pi 



1 + 



6 T 1 3/2 



5 T\ 



(35) 



(36) 



in good agreement with (1+5) = (T v i r /Ti) 3 ^ 2 from equation (30) considering 
that the factor 1/2 in the definition of T v i r in equation ( |35[ ) was somewhat 
arbitrary. In reality, the gas evolution might not be completely adiabatic 
during the collapse, because of the above-mentioned cooling processes. [] We 

4 Our derivation also neglected entropy generation due to the thermalization of bulk 
kinetic energy. When an object virializes, the infall kinetic energy of the gas is thermalized 
in a virialization shock. Thus some entropy is generated and the pressure of the gas is 
higher (typically by a factor 1-2) than predicted by the adiabatic compression. In any 
event, this entropy generation would only decrease the above 6/5 by a factor of order 
unity and would not change the results substantially. 
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therefore adopt the following procedure in our simulations: the gas density 
is evolved according to the top-hat solution until T reaches T V i r . At this 
point, gas pressure is assumed to halt the collapse, and the gas density is 
held constant for the rest of the run. If the gas overdensity reaches the virial 
value 187T 2 before T reaches T vir , then the density is held constant at this 
value, and the temperature is raised to T v i r (by assumed shocks). 

What happens now, after z V i r ? If the gas is going to be able to collapse 
further and eventually form something like population III stars, the baryons 
must now be able to dissipate energy rapidly through cooling. If this is 
the case, the gas cloud may get dense enough to become self-gravitating, 
which adds further instability to the system and may eventually lead to the 
formation of an extremely nonlinear object like a galaxy. The key question 
is thus how fast the gas in the lump can cool after z v { r . This is the topic of 
the next section. 
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5 NUMERICAL RESULTS 



After a lump virializes, one of two things will happen to it: 

• Enough H2 is produced that it will enter a phase of runaway cooling 
and collapse. 

• Cooling will be so slow that it will remain pressure supported for a 
Hubble time. 

In the former case, we will say that the lump collapses, in the latter case 
that it fails to collapse. If it fails, it will not produce any luminous objects 
that can reheat the IGM, but merely remain as an object resembling a 
small Lyman alpha cloud. Whether a lump succeeds or fails to collapse 
of course depends on cosmological parameters such as h, 0, and f2fc. First 
and foremost, it depends strongly on the parameters M and z v i r . In this 
section, we first give an operational definition of what we mean by collapse, 
and then evolve a large number of lumps numerically to see for which parts 
of parameter space they manage to collapse, summarized in figures I and |. 

The results of two sample runs are shown in figures ^ and [|. Both 
have T v { r = 1000K and the standard CDM parameters Q = 1, fij, = 0.06, 
h = 0.5. In Figure || (with z v { r = 100), collapse succeeds by our criterion 
below. To the left, we see how recombination reduces the ionization fraction 
x sharply at z ~ 10 3 . This weakens the Compton coupling to the CMB, 
and at z ~ 400, the gas temperature begins dipping slightly below the CMB 
temperature (straight line). At z ~ 300, a minute fraction of molecular 
hydrogen is formed via the channel before this reaction freezes out. At 
z ~ 100, density and temperature rise to their virial values. This causes a 
surge in the production of H2 via the H~ channel, producing a molecular 
abundance close to 10 -3 , and this in turn causes rapid cooling. From this 
point on, the curves in the figure are of course irrelevant, as the density 
will rise, causing even more rapid cooling and a density profile that must 
ultimately be modeled with a full 3D hydrodynamics simulation. 

The evolution of a less successful lump is shown in Figure f|, with z v { r = 
10. Here even the molecules produced by the //"-channel around z ~ 100 
are too few to cause significant cooling. The molecules produced in the 
third wave wave of formation, at z ~ z v i r , are unable to cool the cloud 
substantially simply because the density (and thus the cooling rate) has 
become too low. Figure ||]. 
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5.1 The collapse criterion 

We now give our operational definition of failure to collapse. After the lump 
has virialized, we keep the density constant at p v ; r and continue to integrate 
the equations for the time evolution of temperature, ionization fraction and 
molecule abundance. Loosely speaking, we consider the cloud a failure if 
its temperature has not dropped substantially within a Hubble time, which 
roughly corresponds to the redshift dropping by a factor 2 2 / 3 . We define 
failure to mean that 

T( V z) > r,T(z), (37) 

and choose rj = 0.75. We do not want to choose r\ too small, since then even 
clouds that merely "loiter" for a while and suddenly cool at a substantially 
lower redshift (when molecule formation suddenly becomes effective) will be 
counted as successful. 

It should be noted that Compton cooling alone is useless for making early 
structures. If it is able to cool the cloud substantially, the resulting contrac- 
tion will drive up the recombination rate (since the CMB temperature is 
<< 10 4 K), virtually all free electrons will disappear, and Compton cooling 
will cease. Thus Compton cooling is self-destructive. Molecular cooling does 
of course not suffer from this problem once the H2 has been produced, and 
can make runaway contraction proceed over many orders of magnitude. The 
same goes for hydrogen line cooling: although it requires free electrons, the 
latter will be produced collisionally at the high temperatures ~ 10 4 K where 
line cooling is effective. 

To ensure that our minimum mass is that above which runaway collapse 
(and thus possible formation of luminous objects) can occur, we thus ignore 
Compton cooling when z < z v i r . 

5.2 The "shooting" scheme 

For each virialization redshift z v i r , there will clearly be some critical tem- 
perature T v i r such that clouds with T > T„j r will collapse and clouds with 
T < T v i r will fail. We find this critical value by a "shooting" scheme: we run 
the code for a very high and a very low virial temperature, then again for 
the average of the two temperatures, then use the interval halving method 
to recursively home on on the critical value T v i r . This is quite feasible nu- 
merically, since each individual evolution run takes merely a few seconds on 
a workstation. 
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6 DISCUSSION 



Earlier work on H2 formation in the early universe has focussed on photodis- 
sociation and subsequent suppression of H2 cooling near the first structures 
to form, which are a likely source of ionizing photons at high redshift (Silk 
1985; Efstathiou 1992). Conversely, Haiman, Rees & Loeb (1996) find that 
at low redshift the ionizing background radiation field from the first collaps- 
ing systems may actually stimulate H2 formation and cooling in primordial 
clouds. We have shown that even without ionizing radiation, enough elec- 
trons survive from the recombination epoch, even in overdense collapsing 
regions, for H2 formation and cooling to be significant. Indeed, we have 
found that in the context of a standard CDM model, H2 formation triggers 
cooling in virialized clouds and allows early formation of low mass objects. 

Typical initial conditions for the first bound objects to form (from 3- 
a peaks) at z ~ 30 are found to be fn 2 ~ 10~ 3 , n# ~ 10 2 cm~ 3 . Clouds 
of baryonic mass ~ 10 5 M Q can be virialized at this redshift, with ensuing 
runaway H2 cooling. The abundance of such objects is readily estimated 
by combining the Press-Schechter formalism with the accurate derivation 
of the small-scale transfer function given by Hu &; Sugiyama (1995). The 
impact of these first objects depends strongly on unknown quantities such as 
their star formation initial mass function (IMF). The most important issue 
is whether they were able to emit enough ionizing radiation to reionize the 
universe or not. Below we will argue that they might have left observable 
imprints in both cases. 

6.1 If UV-emission is substantial... 

Let us first consider the former case, where UV-emission is substantial. Our 
population of condensed baryonic clouds could either undergo star formation 
or form massive black holes. If the former fate awaits these clouds, it is 
plausible to believe, by analogy with our knowledge of the most metal-poor 
Galactic stars, that a wide range of stellar masses is generated. In either 
substantial production of ionizing photons is likely. In the former 
case, heavy elements will also be synthesized. This would give a possible 
source for the heavy elements found at z = 2 — 4 in Lyman alpha forest 
clouds, the most primitive objects in the universe, that amount to ~ 0.3% 
of the solar abundances. The IGM will be reheated at z > 10, thereby 
suppressing the formation of dwarf galaxies until a much later epoch, as 
argued by Blanchard et al. (1992). The low luminosity tail of the luminosity 
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function of faint blue galaxies is indeed inferred to steepen with lookback 
time, as interpreted in models of faint galaxy number counts (Treyer and 
Silk 1994), consistent with recent (z ~ 1) formation. 

In addition, optical depths of at least a few percent (Tegmark and Silk 
1995) to electron scattering in the IGM are inevitable if reionization occurs 
when the first generation of objects condenses. This would lead to notewor- 
thy implications for satellite proposals to measure the CMB anisotropy C[s 
to a precision of a percent or so. Scattering at this level would reduce the 
height of the acoustic peaks, which in the absence of early reionization are 
primarily sensitive to the baryon density. 

6.2 ...and if it is not 

Let us now consider the latter case, where the initial UV-emission is negli- 
gible. Even if star-formation is successful, there are at least three possible 
things that could prevent substantial UV-emission: 

1. The IMF could be so steep that almost no OB-stars are formed. 

2. The bulk of the UV-radiation could be absorbed locally, so that most 
of the radiation leaving the cloud is degraded below the Lyman limit. 

3. Since the clump would be quite loosely bound, with a virial tempera- 
ture <C 10 4 K, the first few massive stars might photoionize the entire 
cloud, blow out the gas and thus prevent the bulk of the baryons from 
forming stars. 

If any of these caveats apply, then a much larger fraction (i.e., not just 
3 — <j peaks) of the baryons would have time to form stars before global 
reionization finally raised the Jeans mass to above 10 4 K and terminated 
this production of small objects. This turn-off might not occur until z ~ 
5 — 10, which could leave as much as 50% of the baryons in condensed 
MACHO-like objects. For a low density CDM cosmology with Q ~ 0.3 and 
a nucleosynthesis-favored baryon fraction Qf, ~ 0.06, this would imply that 
about 10% of our Galactic halo would consist of MACHOs. 

The authors would like to thank Martin Haehnelt, Uffe Hellsten, Avi Loeb 
and Ned Wright for useful comments. This work has been partially sup- 
ported by European Union contract CHRX-CT93-0120 and Deutsche Forschungs- 
gemeinschaft grant SFB-375. 
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APPENDIX A 



In this appendix, we provide fits to the CMB photodissociation rates of H~ 
and H ^ '. 

The cross section for photodissociation of H~ of Wishart (1979) is well 
fit by the expression 

( x - 1)3/2 

a » 3.486 x 10~ 16 cm 2 x v ( — , (38) 

where x = hv/0.74eV. 

To accurately compute the photodissociation rate of , one would have 
to include photodissociation from all its excited states. However, due to lack 
of reliable molecular data, we only use the rates for photodissociation from 
the ground state computed by Stancil (1994), which we find to be well fit 
by the expression 

a « 7.401 x 10- 18 cm 2 x kt* 2 -0-0302^-0.0158^ (3g) 

where x = 2.762 ln(/tz//11.05eV). The cross section vanishes below the bind- 
ing energy hv = 2.65 eV. (Neglecting dissociation from excited states will 
lead to a slight overestimate of the ^-production though the channel.) 
To obtain the desired dissociation rates k, we simply integrate the above 
cross sections against a Planck spectrum: 

Stt too yio{y)dv 

C 2 Jo eW^-l' ^ 

and fit the numerical results by the simple expressions given in Table 1. k 4 
is accurate to within 10% for the redshift range 40 < z < 2000, and k-j is 
correct to within 50% for 150 < z < 1500. 
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Figure 1: Molecular fraction needed and molecular fraction produced. 
The solid, short-dashed and long-dashed lines correspond to lumps virializ- 
ing at z v i r = 100, 50 and 25, respectively. Only clouds above the downward 
sloping lines (outside the shaded region for z v i r = 100) can cool in a Hubble 
time. The upward-sloping lines show the molecular fraction produced in a 
local Hubble time, so the minimum temperature needed for collapse is that 
where the pair of curves cross (solid dots — lower z v { r require higher virial 
temperature). Electron depletion is the limiting factor above the thin dot- 
ted line, so we see that for z > 50, the results are rather independent of the 
initial ionization fraction. 
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Figure 2: Model for density evolution. 
Our model for the evolution of the baryon number density n{z) is shown 
for models with three different virialization redshifts z V i r , for the case of 
negligible pressure, n first decreases slower than the background density 
(dashed line) according to linear theory, then increases again as the lump 
collapses and virializes, and finally reaches the virial plateau value of 18-7T 2 
times the background density when z = z v i r . 
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Figure 3: Lump evolution. 
The time-evolution of gas in a lump is shown for z v i r = 100, T v i r = 2000K, 
h = 0.5, VL = 1 and Qj, = 0.06. From top to bottom on the right side, 
the curves show the number density n in units of 10 3 cm -3 , the molecular 
fraction /, the temperature T in units of W 6 K, the CMB temperature in 
the same units and the ionization fraction x. 
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Figure 4: Lump evolution. 
Same as previous figure, except that z v i r = 10. 
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Figure 5: The minimum virial temperature needed to collapse. 
The minimum T v i r for which collapse succeeds is plotted at a function of 
virialization redshift for standard CDM (f2 = 1, flf, = 0.06, h = 0.5). Only 
lumps whose parameters (z v i r ,T v i r ) lie above the shaded area can collapse 
and form luminous objects. The dark-shaded region is that in which no ra- 
diative cooling mechanism whatsoever could help collapse, since T v i r would 
be lower than the CMB temperature. The solid curves show the tempera- 
ture evolution of the uniform IGM and (187r 2 ) 2 / 3 times this value, so above 
the upper line, gas can attain the virial overdensity without problems with 
pressure support. 
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Figure 6: The minimum mass needed to collapse. 
The function M c {z v i r ) is plotted as a function on virialization redshift for 
standard CDM (Q = 1, Qf, = 0.06, h = 0.5). Only lumps whose parameters 
(z v i r ,M) lie above the shaded area can collapse and form luminous objects. 
The dashed straight lines corresponding to T„ r = 10 4 K and T v i r = 10 3 K are 
shown for comparison (dashed). The dark-shaded region is that in which 
no radiative cooling mechanism whatsoever could help collapse, since T v i r 
would be lower than the CMB temperature. The solid line corresponds to 
3 — 0" peaks in standard CDM, normalized to <7g = 0.7, so such objects with 
baryonic mass Oj x 2 x 10 6 M© ~ 10 5 M© can form at z = 30. 
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